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PACS 05 . 45 . -a - Nonlinear dynamics and chaos 

PACS 45 . 05 . +x - General theory of classical mechanics of discrete systems 

PACS 42 . 25 . Dd - Wave propagation in random media 

Abstract. - We suggest that KAM theory could be extended for certain infinite-dimensional sys- 
tems with purely discrete linear spectrum. We provide empirical arguments for the existence of 
square summable infinite-dimensional invariant tori in the random discrete nonlinear Schrodinger 
equation, appearing with a finite probability for a given initial condition with sufficiently small 
norm. Numerical support for the existence of a fat Cantor set of initial conditions generating 
almost-periodic oscillations is obtained by analyzing (i) sets of recurrent trajectories over suc- 
cessively larger time scales, and (ii) finite-time Lyapunov exponents. The norm region where 
such KAM-like tori may exist shrinks to zero when the disorder strength goes to zero and the 
localization length diverges. 



KAM theory [T] predicts that, under generic conditions, 
Hamiltonian systems with a finite number of degrees of 
freedom, N, close enough to an integrable limit exhibit 
quasiperiodic trajectories which are dense on invariant N- 
dimensional tori in a 2A-dimensional phase space. The 
N fundamental frequencies of those trajectories depend 

' on the torus. When some integer combination of the fre- 
quencies vanishes, resonance is obtained, and tori which 

' are resonant or almost resonant generally break up into 
chaotic trajectories when the Hamiltonian is perturbed 
from the integrable limit. Since there are infinitely many 
possible resonances, infinitely many gaps appear densely 
in phase space. As the volume of the gaps drops exponen- 
tially with the order of the corresponding resonance, the 
persisting tori form a fat Cantor set (i.e., of nonvanishing 
Lebesgue measure), which goes to full measure at the in- 
tegrable limit. There is no general extension of the KAM 
theory for infinite systems (N = oo), except for some spe- 
cial models [2]|3]. However, it is generally believed that 
most KAM tori disappear when the dimension of the dy- 
namical system is infinite (see, e.g., [1]). 

Simple empirical arguments confirm that KAM tori, 
which are spatially localized (square summable, Zg) can- 
not survive when the spectrum of the linearized system is 
absolutely continuous in some frequency interval. This sit- 
uation occurs, e.g., for spatially periodic arrays of coupled 



anharmonic oscillators. In such models, any hypothetical 
quasiperiodic solution with at least two incommensurate 
frequencies would generate harmonics densely on the real 
axis, overlapping with the interval of the absolutely con- 
tinuous part. Thus, these harmonics would radiate energy 
towards infinity so that the localized energy could not be 
conserved [5]. (By contrast, simple periodic solutions may 
remain localized, forming intrinsic localized modes ["dis- 
crete breathers"] which may be dynamical attractors for 
some initial conditions [6 ].) However, these arguments do 
not hold when the linear spectrum is purely discrete, and it 
is known, e.g., that spatially localized periodic solutions 
with frequencies inside the linear spectrum exist generi- 
cally in systems with linear Anderson localization [TJ[8] . 

We consider here the one-dimensional random discrete 
nonlinear Schrodinger (DNLS) system (see, e.g., [9HT%]). 



ii>n = (e n + X\^n\ 2 )^n ~ C(lp n+ i + 1p n -l), 



(1) 



with integrable limits at \ = (linear limit) and C = 
(anticontinuous limit). The random onsite energies e n 
are uniformly distributed in the interval [— W/2 : W/2]. 
(See [H] for a direct experimental realization of (JlJ with 
optical waveguide arrays.) Then the linear spectrum 
(x = 0) is discrete and the corresponding eigenstates 
are lg (exponentially localized). Bourgain and Wang [16] 
have proven that close enough to the linear limit, eq. ([TJ 
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exhibits quasiperiodic solutions corresponding to finite- 
dimensional tori in phase space. Since in finite systems, 
invariant tori with dimension P < N generally have zero 
measure, we expect similarly the solutions found in |16) to 
have zero probability to occur in an infinite system. In this 
work, we provide empirical arguments, supported by nu- 
merical calculations, suggesting that eq. ([IJ may also sus- 
tain infinite dimensional invariant tori of almost periodic^] 
trajectories, which are 1% and could be found with nonva- 
nishing probability for an initial condition at small enough 
norm. Frohlich, Spencer and Wayne proved the existence 
of such tori [2J, but only in a special family of modified 
random DNLS-type systems not including eq. ([T]). 

Our empirical arguments can be summarized as fol- 
lows: We consider an arbitrary (g small-amplitude ini- 
tial condition, corresponding to a distribution of excited 
linear modes coupled by weak nonlinearities. Estimating 
the probability that these couplings induce resonances be- 
tween the linear modes, we find this probability to vanish 
in the limit of small norm. Since our argument explic- 
itly uses the fact that the linear spectrum is discrete with 
exponentially localized eigenstates, it does not hold for a 
system with an absolutely continuous part in its spectrum. 
Moreover, we predict that the existence domain of KAM- 
like tori shrinks to zero in the limit of weak disorder, or 
equivalently long localization length. 

For the linearized eq. (JTJ) , the discrete (countable) set of 
eigenvalues ui p is associated with a basis of real I2 eigen- 
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panding {ipn} in this basis, ip n (t) = 



n — J 



Ex- 
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and the norm square 



which is the second conserved quantity of eq. ((TJ (in ad- 
dition to the Hamiltonian) . Then, we obtain for the new 
complex coordinates \i p (cf, e.g., [91 IT51IT3] ) . 



ifi p = u p fi p + ^2 C p,p' CO/V (*)> 



(2) 



p' 



where C p>p ,(t) = xErfCjCWWV,!,! 1 are defined 

via the overlap sums V P , v ', q , q ' = Xl n W W 4>n 4>n ■ 
The coefficients C PtP /(t) are real and depend on time 
through the coordinates themselves. 

In the limit of small amplitude \fJ- p \, cubic terms in 
eq. ([2j are higher order and may be neglected during some 
time. The linear behavior is n P (t) ~ (j v (0)e~ ta>t ' t . Then 
during that interval of time, coefficients C PtP i(t) are almost 
periodic functions of time (i.e., with a countable set of pe- 
riods), C PtP >(t) « xE M ^ P y, M 'K(0)^(0)e ,(u ' 



,)t 



1 An almost periodic function /(£) can be written as an absolutely 
convergent series, /(£) = JZ n / E e'""', where the set of frequencies 
u) n is countable and ^ |/ n | < +00. |17l 



Fromcq. ©, || Mp | 2 = 2 J2 p ,^ p Cp,p,(t)Im(^ pf i p ,), defin- 
ing the norm current J 2 
son modes p 7^ p' as 
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between two different Ander- 



J p >^ p = 2C PtP >{t)lm(n*[i p >) 
2xE [^.^K(0)l ' IM<>)I eKK-^-K-v)) 
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x|/Xp(0)| • |/v(0)l sin((o;p - u p >)t - [a p - a p >)), (3) 



where a p is the initial phase of jU p (0) = |^ p (0)|e _lap . This 
current is almost periodic in time, and its oscillations 
should be small in order that |Mp(i)| 2 remains approxi- 
mately constant, so that /%,(£) « /U p (0)e -I " p * remains valid 
for all times. Since time integration of the current yields 
denominators ui q — u> q i ± (u) p — uj p > ) which may be small, 
this condition requires that there should be no strong reso- 
nances between any pairs of sites p ^ p' '. We consider res- 
onances to be weak enough when the numerator is smaller 
than the denominator for all terms, i.e., 

|w g -uy±(w p -uy)| £ «Ix^, 9 ,?'Hm<z(0)H/v(0)I ( 4 ) 

for all q and q' , where k is of order 1 (Chirikov criterion). 
Note also that ± can be dropped, since the same condition 
is obtained if q and q' are exchanged. 

We assume now in order to fix the ideas that io p are ran- 
dom numbers distributed in some interval with a smooth 
probability law with maximum density Pq. Each reso- 
nance p,p',q,q' has thus a probability with the upper 
bound 2P n\xV p y >M /| • |/x g (0)| • |/v(0)| to occur. The 
probability Pr that there is at least one resonance in the 
system for this initial condition, is bounded by the sum of 
these probability bounds divided by 2 (resonance p, p' , q, q' 
is equivalent to p' ,p, q',q). 

Pfl<p<HxlEK(o)K 9 'l/v(o)|, 

9,9' 

where A q y = E p ^ p ' m>,p',q,g'\- Then we obtain 

p fl <Po«lxH|A||.||W0)}|| 2 , 

with A = {A qq i} and the norm 

IIAII ll A ' X IH 

All=aup " .. ■■ '" . 
x ll-X-ll 

||A|| may be equivalently defined as the smallest upper 
bound of the eigenspectrum of A. 

The probability to have at least one resonance is thus di- 
rectly related to the norm square | |{/J-g(0)}| | 2 of the initial 
condition (this was also found numerically in |13j). The 
probability P/v = 1 — Pr to have no resonance thus has a 
tower bound, Pjv > 1 - Po«|xl ' ||A|| • ||{M g (0)}|| 2 , so that 
when the norm of the initial condition is small enough, 



IIWo)}|| 2 < 
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we obtain P/v > 0. Thus, if we can prove that the norm 
||A|| is not infinite, the probability to have no resonance 
will be nonvanishing, and go to 1 when the norm of the 
initial condition goes to zero. 

An upper bound for ||A|| can be obtained from 



|A||<sup]T|A 



<?,<? 



<su P J2 El^M^H^H^I 

q pjtp',q' n 

<sup^(^|^)|) 3 |^)|. 

a 

^ n p 

If the eigenstates are exponentially localized, then 

Yj P \4>n | < +00 and J2n l^« I < +°°- ^ we assume 
these are bounded for all n or q by the same constant S, 

then 1 1 A|| < S 4 < +oo@ Note that since £ n |4 9) | 2 = 1, 
if the localization length increases and diverges, then 
S > J2 n \vn I — > 00. To obtain a reasonable estimate of 
||A|| we assume an exponential bound for all eigenstates, 



<K 



1-A 2 
1 + A 2 ' 



\\n-P\ 



where K is some constant and A = e -1 ^, where £ is the 
localization length. Then S* 4 = K 4 (\\2\2 n~rn - When 
the localization length diverges at weak disorder we find 



l|A|| < 16K 4 £, 2 . 

Consequently this upper bound for the norm diverges, sug- 
gesting that ||A|| might also diverge in the same way. 

We did not consider the probabilities of resonances at 
higher orders 6,8,... which are cumbersome to calculate. 
The correction on the bound of Pr would be higher or- 
der in 1 1 {^^(O)} 1 1 2 and thus could be neglected in the 
limit of small norm. We conjecture that at each order 
2p, these probabilities can also be bounded by convergent 
series multiplied with ||{^g(0)}|| 2 ( p_1 ). The probability of 
higher order resonance is expected to decay exponentially 
with the order. Thus when the norm of the initial condi- 
tion is not too large, we would expect that the probability 
of no resonance at any order is non-vanishing and still 
bounded from below, going to 1 as the norm goes to zero. 

Note that for finite systems with size N, the linear spec- 
trum is always discrete and the series for bounding Pr 
become finite sums, implying ||A|| < +00. Then, we 



2 If /i. g (0) is not arbitrarily chosen, we may get a better upper 
bound for the existence of KAM tori. For example, if li q (0) = 
$q,q o \\{fiq(0)}\\ is initially localized at a single Anderson mode 
q' = go [I3|, we have P R < P K.\x\(sup q , \A q ^ q ,\)- |{> g (0)}|| 2 where 
sup / \A„i „i\ < ||A||. Thus, initial wave packets which are close to 
single Anderson modes survive much better as almost periodic solu- 
tions than those which are arbitrarily spread (in Anderson space) at 
the same norm. This effect is especially important when the local- 
ization length is large since then sup / |j4 ? / i9 /| << ||A||. The same is 
true if the initial wavepacket is split into several packets with smaller 
norm which are far apart at the scale of the localization length. 



know that the conclusion of our empirical argument is 
consistent with KAM theory, predicting the existence of 
A— dimensional invariant tori of quasiperiodic solutions 
at small enough amplitude (or equivalently, for finite sys- 
tems, small enough norm) with a probability going to 1 at 
zero amplitude. Our conjecture is that this argument also 
holds for infinite systems provided ||A|| < +00. This sit- 
uation occurs when the linear spectrum is purely discrete 
with exponentially localized eigenstates, but is not fulfilled 
when it contains an absolutely continuous part. Then we 
would conclude, that the norm region for initial conditions 
where KAM tori may exist shrinks to zero when the lo- 
calization length diverges, approaching the limit without 
disorder where the linear spectrum is absolutely continu- 
ous and no 1% almost periodic exact solution could survive 
due to radiaton. 

Attempting to distinguish numerically KAM tori among 
other solutions in the random DNLS equation, we first 
consider that trajectories of KAM tori are almost periodic 
in time, and use the Harald Bohr theorem: 

If a function f(t) is almost periodic, then for any e > 0, 
there exists a relatively dense set of translations r such 
that \f(t) — f(t + r) I < e for all t €] — 00, +00 [. In other 
words, there exists a diverging subsequence r n fulfilling this 
property with r n < T n +i and r n +i — r„ bounded. 

Thus, for KAM-like tori, recurrences should be observed 
numerically in all quantities which depend on time such 
as local coordinates at arbitrary sites, momenta, partici- 
pation number, etc. Typically, if recurrences appear for 
one of these quantities, they are found also for any other. 
Some problems with this method are obviously that 

(i) the Harald Bohr theorem can be checked only for 
finite times r £ [0,T] and t € [0,T], where T is the time 
of integration, and for finite-size systems; 

(ii) then e cannot be chosen too small for avoiding that 
recurrences become too rare, and thus that the corre- 
sponding pseudoperiod exceeds the integration time; 

(iii) the system size should be sufficient in order that the 
amplitude is practically zero at the edge, so that boundary 
effects can be neglected during the integration time T. 

(iv) the integration accuracy should be good enough, 
for avoiding numerical drift from KAM-like tori to neigh- 
boring chaotic trajectories. In practice, the relative er- 
ror in the conserved quantities are kept at the order 10 -6 
or smaller in all simulations, with consistency checks for 
smaller systems reaching accuracies 10 ~ 8 — 10~ 10 . 

If there are KAM tori persisting over infinite time, we 
should expect that their structure as some model param- 
eter varies is a fat Cantor set with infinitely many gaps 
due to resonances at all orders, but nonvanishing Lebesgue 
measure (or nonvanishing probability). Thus, one should 
find that the probability (in the space of initial conditions) 
that recurrence persists over a time T does not shrink to 
zero as T becomes very large. However, high order reso- 
nances (corresponding to small gaps) should manifest only 
after very long integration time. These trajectories may 
look almost periodic and exhibit recurrence over relatively 
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disorder strength W 
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Fig. 1: (a) Last observed time for recurrence (e/N = 0.02) 
in I'i/'nol 2 versus norm in a particular disorder realization with 
W = 20. (b) Sets of T-recurrent trajectories at T = 10 4 for 
various disorder strengths W . e/N varied from 0.1 (W — 6) to 
0.005 (W > 14). (c) Fraction of the trajectories T-recurrent at 
T — 10 3 (excluding the high-norm self-trapped regime) which 
remain T-recurrent also at longer times, for various disorder 
strengths. At final integration times, e/N is, respectively, 0.03 
(W = 6), 0.003 (W = 12), 0.001 (W = 30), and 0.0005 (W = 
20). In all figures, the disorder realization is the same, the 
initial condition ip n (t = 0) = VNS n ,n with e„ « 0.465291F, 
C — 1, x — ~ 1] an d system size N = 500. 



long times, before they blow up as chaotic trajectories. 
(Such trajectories were numerically identified in |13] as 
belonging to a "regime I" of rather small norm and/or 
strong disorder.) Generally, all trajectories which remain 
recurrent after a time T will be termed T-recurrent. 

An illustration is given in fig.[TJa), showing the last ob- 
served recurrence time as a function of norm for a single- 
site initial condition in a particular disorder realization for 
rather strong disorder W = 20 (we chose the initial-site 
energy e na rather close to the upper band edge, in order 
that an increasing negative nonlinearity will "scan" most 



possible resonances inside the band). Any horizontal in- 
tersection of this graph at a given time T yields the set 
of initial norm generating T— recurrent trajectories. As 
can be seen, many trajectories remain recurrent for times 
larger than 10 7 . In addition to the regime of small norm 
(here N < 3.6) expected from our argument above, there 
is also a regime of recurrent states for N > 19, as well as 
a small interval around N ~ 15.6. The recurrent trajecto- 
ries in the large-norm regime are related to the fact that, 
above some threshold norm, an increasingly larger part 
of the norm will self-trap around the initial site no [T2"] . 
The total norm available for the rest of the lattice will 
then actually decrease, and essentially the same argument 
as above could be used in support of a KAM-like regime. 
Note also that the simultaneous limit of strong disorder 
and large norm, W — > oo,A/" — s- oo with X-A/"/W finite, 
is equivalent (rescaling time) to the anticontinuous (un- 
coupled) limit C = which is also integrable, and where 
the linear localization length is vanishing. In this limit, 
for any initial condition, each anharmonic oscillator n ex- 
hibits a periodic motion with frequency u>' n ^ e n different 
from the corresponding onsite linear frequency. We may 
thus expect a KAM-like regime also close to this limit. 

Thus, though numerics cannot provide a rigorous proof, 
the most plausible interpretation is that that there is 
an underlying Cantor set of initial conditions generat- 
ing KAM tori, persisting over infinite time and infinitely 
large systems. Note also that the gap structure of Fig. 
QJa) is reminiscent of the "stickiness" phenomenon in low- 
dimensional systems, where many initial conditions close 
to (but outside) the Cantor set of KAM tori remain close 
for long times before they finally escape (compare, e.g., 
with Figs. 7 and 10 in [T5] for the standard map). 

The variation of the sets of T-recurrent trajectories 
with the disorder strength is illustrated by fig. QIb), for 
a rather modest time T = 10 4 . As predicted, there is 
always a small-norm regime where most trajectories are 
T-recurrent, the size of which grows with increasing disor- 
der strength. There is also always a high-norm T-recurrent 
regime with lower boundary increasing with increasing dis- 
order, since the norm necessary for efficient self-trapping 
increases linearly with W for a single-site initial condition 
[12"] . In-between these two regimes, for larger W there 
are also several intermediate regimes of T-recurrence, sep- 
arated by gaps with non-recurrent (chaotic) trajectories. 
For smaller W the relative sizes of these gaps grow, and 
they merge into one single main gap of trajectories which 
typically show a chaotic time-evolution and spread sub- 
diffusively [9"r [m[l"3] (for long but possibly finite times). 

Similar pictures as fig. Hfb) can be obtained for larger 
times, although obtaining a good resolution with sufficient 
numerical accuracy to clearly identify recurrences makes it 
very time-consuming for times larger than ~ 10 6 . As a rule 
of thumb, to determine persistent recurrences e is divided 
by two for each order of magnitude in time. In fig. [He), 
we give an example showing how the measure of the set 
of T-recurrent initial conditions varies with T, for differ- 
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ent strengths W of the same disorder realization. In order 
to obtain a finite set, we here exclude recurrent trajecto- 
ries belonging to the high-norm (self-trapped) regime, and 
moreover for comparison we normalize the sets by dividing 
with the number of T-recurrent trajectories at T = 10 3 . 
The data of fig. [He) suggest the existence of an asymptotic 
set with a nonvanishing measure at infinite time. However, 
for stronger disorder, it is clear from fig. [TJc) that to get 
a clear picture of the asymptotic measure of this set, con- 
siderably longer integration times would be needed. This 
is due to the fact that the fraction of "sticky" trajectories, 
which are T-recurrent over very large but finite times, in- 
creases with the disorder strength (cf. fig. HJa)). 

Thus, the numerical study of T-recurrent trajectories 
for finite times can give us only an indication about the 
true nature of the KAM-like trajectories, and in particular 
the presence of long-time "sticky" trajectories makes it ex- 
tremely difficult, e.g., to resolve the Cantor-set structure 
of resonances in the low-norm regime within a reasonable 
amount of computer time. We therefore now turn to dis- 
cusss another technique to numerically distinguish KAM 
tori, in terms of the tangent map and the corresponding 
finite-time Lyapunov exponents. Small perturbations of 
eq. ([U yield the Hill equation 

il'ln = (e„ + 2x\lpn\ 2 )Vn + X^Vn - C(T] n+1 + J/ n _l). (5) 

If eq. (p} possesses an almost periodic solution with the 
discrete set of frequencies u>i,W2j ■•• ) w p , ..., corresponding 
to KAM tori with full dimension, 

ip n = F n (u!it+ai,u 2 t+a2, ... ,u p t+<Xp, ...;ui,u>2, —,u p , —)>| 
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where F n {x\,X2,...x p ,...;{jJi,U)2,...,U) v ,...) is 27r-periodic 
with respect to xi,X2,-..,x p ,..., then T] n {t) — -g^r- is an 
almost periodic solution of eq. ([5]) . One gets a complete 
basis of solutions of eq. ([5]) by adding also the solutions 
Vn(t) = %£ = t%* + §§J, where also |^ is almost peri- 
odic. Consequently, if the solution of eq. ([1]) corresponds 
to a KAM torus, the general solution of eq. ([5]) grows lin- 
early as a function of time. Note also that for any solution 
ip n to eq. (HJ) , there is a trivial solution r\ n = iijj n to eq. (JS|) 
corresponding to a global phase rotation. Numerically, we 
remove this component by subtracting the projection of 
r) n (t), obtained by integrating a randomized initial condi- 
tion rj n (0), on this vector. 

Thus, considering the total norm of the perturbation 
divided by t, j||j7(*)|| where ||7?|| = \ZY, n \Vn\ 2 , this quan- 
tity will exhibit bounded oscillations for all times if ip n 
corrsponds to a KAM torus, and diverge exponentially 
with a positive Lyapunov exponent for any chaotic trajec- 
tory. A numerical illustration is given in fig.[2ja), showing 
a narrow resonance gap in the low-amplitude KAM-like 
regime around N = 0.600 for disorder strength W = 12. 

From the so obtained numerical solutions r] n (t), we 
calculate finite-time Lyapunov exponents as A(ijvf) = 



Fig. 2: (a) Total norm of perturbation T} n (t) divided by time, 
for three solutions corresponding to single-site initial conditions 
tp n with slightly different N ~ 0.6. (b) Finite-time Lyapunov 
exponents for 8 solutions in the same regime as in (a). At time 
10 6 , the upper curves in (b) correspond, from top to bottom, 
to M = 0.604,0.600,0.602,0.601, while the lower curves for 
TV = 0.595, 0.599, 0.603, 0.605 all follow very closely a curve ~ 
logt/t, as expected for KAM tori. Disorder strength W — 12, 
other parameters and disorder realization same as in fig. [T] 
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r, where to = 0, and t m are chosen 



to correspond to a Poincare section defined by |^>n (£)| 2 
having a local maximum at each t = t m . Thus, for a re- 
current trajectory, the optimal recurrence times t^ form 
a subset of the t m . Moreover, for an almost-periodic tra- 
jectory with a time-linear growth of \\rj\\, A(t) should de- 
crease asymptotically to zero as A(t) ~ log(t)/t for large 
times. As can be seen from the example in fig. [H[b), the 
behaviour of A(i) in the neighborhood of sharp resonances 
is very sensitive to small parameter variations. In this ex- 
ample, the trajectories are seen to be (possibly weakly) 
chaotic in the interval 0.600 < J\f < 0.602 and in an even 
narrower interval around Af ~ 0.604, while apparently 
almost-periodic trajectories (with no visible deviations 
from the asymptotic behaviour A(i) ~ log(i)/i for times 
larger than 10 6 ) are seen e.g. for N = 0.595, 0.599, 0.603 
and 0.605. Thus, this again supports the existence of a 
finite-measure Cantor set of almost-periodic KAM tori. 

Comparing the numerical results from the analysis of 
recurrences and tangent map, they are consistent in the 
sense that when recurrences are lost, there is a clear devia- 
tion in 1 1 r\ 1 1 from time-linear growth, and in most cases the 
growth is exponential with well-defined non-zero A(i) (at 
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Fig. 3: Ensemble-averaged fractions of trajectories T-recurrent 
at T — 2 ■ 10 5 versus chain size, for random single-site initial 
conditions with < N < 20 and disorder strength W = 12. 
100 different realizations were used for each size, and dashed 
lines show the standard deviations of the distributions. Trajec- 
tories were considered T-recurrent if A(T) < 10 -4 (cf. fig. |2J. 



with nonzero probability: (i) those generating spatially lo- 
calized (non-spreading), almost periodic solutions (KAM 
tori); and (ii) wave-packets which are initially chaotic and 
spreading [5 HTT][T5] . Our results should stimulate further 
attempts towards more rigorous treatments, as well as 
more detailed numerical studies of the role of such tori in 
various situations. Ongoing work suggests, e.g., that also 
large-norm wavepackets extended on many sites could gen- 
erate KAM tori with high probability if the norm density 
is small enough. Thus, the spreading of chaotic wavepack- 
ets may after long times drastically slow down, if they 
become sticky to such tori. The KAM-like trajectories 
may also be experimentally observable, e.g., in disordered 
waveguide arrays 15 , or in disordered bosonic systems 



least for long times). However, generally (and in partic- 
ular for strong disorder with many "sticky" trajectories), 
the tangent-map criterion is considerably more sensitive, 
and may signal a chaotic behaviour several orders of mag- 
nitude in time before recurrence is finally lost. 

Generally, we find the KAM-like trajectories to be ex- 
ponentially localized, decaying essentially with the largest 
linear localization length £ (e.g, £ ma x ~ 1.1 for W = 12 
as in fig. (2). Thus we may expect that, for stronger disor- 
der, the KAM-like trajectories observed for large systems 
should be only slight perturbations of true KAM tori ex- 
isting for short chains. This was also confirmed numeri- 
cally. An illustration is given in fig. [3] where, for W = 12 



that e„ also is chosen randomly: if e„ is close to the 
linear band bottom, self-trapping occurs already for small 
values of Af, and therefore most trajectories will be KAM- 
like. By contrast, realizations with e„ close to the band 
top as in figs. Q3[2] should yield the largest amount of reso- 
nances. Indeed, for the particular realization used for the 
large-size simulations in figs. [I|[2] we obtain for W — 12 
a total fraction of 38% T-recurrent trajectories in the in- 
terval < N < 20 after T = 10 5 (decreasing to 37% at 
T = 10 6 ), essentially coinciding with the lower bounds in 
fig- El (i.e., this realization does not behave exceptionally). 
In summary, our empirical and numerical arguments 
suggest that in the infinite random DNLS model ((T|), 
there are two kinds of initial wavepackets both occuring 
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